Interplay of charge, spin, orbital and lattice correlations in 
colossal magnetoresistance manganites 
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We derive a realistic microscopic model for doped colossal magnetoresistance manganites, which 
includes the dynamics of charge, spin, orbital and lattice degrees of freedom on a quantum me- 
chanical level. The model respects the SU(2) spin symmetry and the full multiplet structure of the 
manganese ions within the cubic lattice. Concentrating on the hole doped domain (0 < x < 0.5) 
we study the influence of the electron-lattice interaction on spin and orbital correlations by means 
of exact diagonalisation techniques. We find that the lattice can cause a considerable suppression 
of the coupling between spin and orbital degrees of freedom and show how changes in the magnetic 
correlations are reflected in dynamic phonon correlations. In addition, our calculation gives detailed 
insights into orbital correlations and demonstrates the possibility of complex orbital states. 

PACS numbers: 71.10.-w, 71.38.-k, 75.47.Gk, 71.70.Ej 



I. INTRODUCTION 



The obset^ation of the colossal magnetoresistance ef- 
fect (CMR)EI'qB in doped manganese oxides with per- 
ovskite structure, Ri-^A^MnOs (R = rare-earth, A = 
alkaline-earth metal), moved thes&jnaterials into the fo- 
cus of intense research activitjoEHJ. It turned out soon 
that the complex electronic and magnetic properties of 
the manganites depend on a close interplay of almost all 
degrees of freedom known in solid state physics, namely 
itinerant charges, localised spins, orbitals, and lattice vi- 
brations. On the one hand, the strong Coulomb interac- 
tion U and the Hund's rule coupling Jh introduce a spin 
background ^ind affect the charge mobility via double- 
exchangaJ'ElEI. On the other hand, the cubic environment 
of the Mn sites within the perovskite lattice results in a 
crystal field splitting of Mn-d-orbitals into e g and t2 g and 
gives rise to an orbital degeneracy in the ground-state of 
Mn 3+ ions. This orbital degeneracy, in turn, connects 
the electronic system to the lattice, making-ii sensible to 
Jahn- Teller distortion and polaronic effectsEHl. 

There are numerous attempts to describe the electronic 
(i.e., charge, spin, and orbital) interactions of the man- 
ganites on a microscopic level (see, e.g., Refs. |lT| , |l2 13). 
However, most of these models do not reflect the full 
multiplet structure of the Mn ions caused by the cubic 
site symmetry or violate the spin Rotational invariance. 
To our knowledge, Feiner and OlesEil were the first who 
derived a consistent spin-orbital model for the undoped 
compounds. In the following we extend their derivation 
to the case of finite doping and complete the resulting 
electronic model with dynamic electron-lattice interac- 
tions of Jahn- Teller and breathing type. Up to now 
models of similar complexity have only been discussed 
on the basis of rather extensive approximations, where 
mean- field approaches are probably the most popular. 
In particular, the lattice was always treated within the 
adiabatic approximation which is questionable in view of 



the comparable energy scales of spins and phonons. For a 
better understanding of this complex many body system 
we prefer to consider all interactions on an equal foot- 
ing. Using exact diagonalisation techniques we study the 
ground-state properties of the derived model on a four 
site cluster. 

For the hole doped region (0 < x < 0.5) our calcu- 
lations show how the lattice can effectively control the 
spin and orbital correlations and the charge mobility. In 
the undoped compounds the electron-phonon interaction 
suppresses the coupling of spin and orbital degrees of 
freedom and is most effective in determining orbital or- 
der. At doping x = 0.25, where ferromagnetism is sta- 
bilised by the double-exchange mechanism, the coupling 
to the lattice can cause self-trapping of the charge carri- 
ers, which immediately switches the spin order. In turn, 
this change is reflected in dynamic lattice correlations, 
namely the bond length fluctuation. For the half-doped 
manganites we find that the electron-phonon coupling 
enhances the susceptibility for charge ordering which 
is present already due to the strong Coulomb interac- 
tion. In addition, moderate variations in the strength of 
the electron-lattice coupling trigger pronounced changes 
in spin and orbital correlations. This relates to the 
complicated patterns of spin, charge and orbital order- 
ing observed experimentally for different manganites at 
x = 0.5. 



II. MICROSCOPIC MODEL 
A. Electronic interactions 

The electronic and magnetic properties of the mixed- 
valence manganites are governed by the manganese d- 
electrons, which are divided into <2g and e g according to 
the cubic symmetry within the crystal field. Due to the 
large Coulomb and Hund's rule interactions each of the 
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three t 2g and the two e g levels carries at most one electron 
and the spins of several d-electrons are aligned in parallel 
at every site (see Figure [D) . Starting from the undoped 
compounds with only Mr? + ions doping will remove e g - 
electrons. Hence, the local electronic Hilbert space can be 
restricted to the large Hund's rule ionic ground-states of 
manganese ions in a cubic crystal field (see Ref . [l5]) . For 
Mn 3+ (d 4 ) this corresponds to the spin-2 orbital doublet 
5 E [t 3 2 ( 4 A 2 )e], 
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Here the operators v aa create spin a electrons in the 



(a 



?, e) or t 2g (a = £,r),Q orbitals and S 1 " 1 



c lt c ai raises the on-site spin. The overlap of neigh- 
bouring manganese d and oxygen p orbitals allows for a 
tunneling of the manganese electrons between adjacent 
sites. Due to the specific symmetry of the,iavolved or- 
bitals this hopping acquires an anisotropyE3E3, 
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which can be expressed in terms of the spatial rotations 



R x = (ClY , R y = (Cf ) 2 = (Clr 1 , R z = (Cl) 3 
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The transfer amplitudes t and t w of e ff and t 2g elec- 
trons are small compared to the local energies U and 
Jh- Therefore the mobility of the e g electrons depends 



FIG. 2: Second order virtual excitations contributing to the 
model Hamiltonian H e i, Eq. (|s]). The shaded region corre- 
sponds to t 2 terms, the other terms are proportional to t^. n 
denotes different nonnegative prefactors of A (see Eq. 



on the correlations of the spin background formed by the 
t 2g electrons. Their kinetic energy is maximal, if neigh- 
bouring spins are aligned ferromagnetically, which is 
essence of the well known double-exchange interactior 
Second order processes in t and are responsible for nu- 
merous superexchange interactions between the localised 
spins, which inherit the orbital anisotropy of the hopping. 

The electronic part of our model is derived by sec- 
ond order degenerate perturbation theoryO, i.e., we cal- 
culated the matrix elements of H t between the—basis 
states (0) or (|^) and all admissible excited statealj. In 
Figure pTthe spin part of the fifteen different virtual exci- 
tations contributing to the superexchange is summarised 
graphically. We neglect terms which involve three differ- 
ent lattice sites. 

A compact expression for the resulting electronic 
Hamiltonian, 



H e i - Rs(Hf i+s + H.c.) 



(5) 
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FIG. 1: Local electronic structure of Mn-d-electrons in a cu- 
bic environment. Here A c f, Jh, and -Bjt denote the crystal 
field splitting, the Hund's rule coupling, and the Jahn- Teller 
splitting, respectively. 



is obtained by rewriting the orbital and charge degrees 
of freedom of the basis states (|l|) and (||) in terms of new 
Fermi operators (V a and projectors P a , 

|0> = 4|O>, |e> = 4|0>, |a 2 > = 44|0>, 

P = n d (l-n £ ), P" = n e (l-n e ), P a *=n e n g . (6) 

Note that these projectors are related to the common 
pseudo-spin operators 



1 ~ z 
2 a . 



r X /V — 



R x /y(T z ) = i(-a* T V3* x ), (7) 



by the equation Rs(P e / E ) = | ± t 5 . Since the amplitude 
of the on-site spin is different for the basis states ([!]) 
and (Eh, it is convenient to represent the spin degree of 
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freedom by Schwinger bosons and which allow for 
a uniform description. Then, the spin operator is defined 
by 2 S = a)^cr^ v a vl but the amplitude is subject to the 



constraint 2 |S| = aja^ + = 4 — n g n E . Using this 
notation in z-dircction the interaction between nearest- 
neighbour sites is given by, 
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(8) 



The first order in i cgjsesponds to the well known double 
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exchange interactionO'Q'EI, which the authors discussed in 
detail recentlyta. The second order in t and £„■ appears 
to be more involved, since the number of accessible vir- 
tual excitations is rather large (see Figure ||). However, 
in all cases it is basically the product of a Heisenberg- 
type spin interaction and two orbital projectors. The 
energies of the virtual excitations depend, in general, 
on all three Racah parameters^ A, B, and C. Apply- 
ing the convention of Ref. |lj we define U and Jh in 
terms of d A d A ^ d 5 ( 4 A 1 )d 3 and d 4 d 4 ^ d 5 £L4i)d 3 ex- 
citations. Together with the approximatiorO C m 4B 
this yields U = A + 22B and J h = 6B. Typical esti- 
mates for thes||^jaxgies are U ~ 6 ... 9 eV and Jh ~ 



B. Electron-phonon interaction 

Concerning the electron-lattice interaction we concen- 
trate on the local environment of the manganese ions, 
namely the surrounding oxygen octahedra (see Figure ||) . 
At every site two optical phonon modes of E g symmetry, 
qs and q £ , couple to the orbital degree of freedom of the 
e g electrons. In addition, a breathing-mode q ai is sen- 
sitive to the electronic density. To lowest order in the 
elongation q a — tf a + b a (a e {9, e, o-iiX this is modelled 
by the E ® e Jahn- Teller HamiltoniarO 

Hjt =9^2 ~ n i,e)( b l,e + he) 

i 

+ (dl g d h£ + dld he )(bl + b h£ )] (9) 
and a Holstein-typeEl interaction 
Hbr =9^2{n it e+n i e - ^n i e n i e ){b\ ai +b i<a J. (10) 






FIG. 3: Distortions of the oxygen octahedra for the Jahn- 
Teller (qe, q e ) and the breathing-type (q ai ) phonon modes. 



The bosonic operators b i a describe the usual harmonic 
lattice dynamics. To a good approximation we can as- 
sume these phonons to be dispersion-less, 



H p h — u) 



(11) 



In experimentsEE the vibrational modes involving dis- 
tortions of the oxygen octahedra are found to have fre- 
quencies of the order of 50 to 80 meV. In the following we 
assume g — g and u> = u. This approximation reduces 
the number of model parameters and is reasonable, since 
we are mainly interested in qualitative features of the 
electron-lattice interaction. 



III. NUMERICAL RESULTS 

The Hilbert space of the complete microscopic model, 
H = H c i + ifjT + Hbr + Hph , (12) 
is large and grows rapidly with the system size. Howeve 



using a density matrix based optimisation procedure^ 
for the phonon subsystem we are able to retain the full 
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quantum dynamics of the lattice and the electronic sub- 
system in our numerical calculation^ of the ground-state 
properties. In addition we take into account some of the 
discrete symmetries of the model, namely the conserva- 
tion of the S z component of the spin, the particle number 
conservation, and the mirror symmetries orthogonal to 
the x and y axes (dot-dashed in Figure ||). Nevertheless, 
the typical dimension of the eigenvalue problem is of the 
order of 10 6 and its repeated solution during the phonon 
optimisation requires the use of large scale computers. 



J h = 0.7 eV, t = 0.4 eV, t/t =3, to = 70 meV 
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FIG. 4: Spatial symmetries of the considered cluster. 



FIG. 5: Left: Total spin S tot and orbital order (oc (n e - n £ )) 
of the ground-state of the cluster in dependence on (a) the 
electron-phonon coupling g, and (b) the Coulomb interaction 
U . Right: Phase diagram of the electronic model without 
electron-phonon interaction (g = 0). 



A. Undoped manganites 

The undoped manganese oxides (LaMnC>3, PrMnOs) 
usually exhibit A-type anti-ferromagnetic order and 
strong J | ah ] tbjTeller distortion of the ideal perovskite 
structureEl'El. The origin of the observed magnetic or- 



J„ = 0.7eV, t = 0.4eV, t/t =3, o> = 70meV 



der has been subject to 
band structure calculations 



;s. While different 
emphasise the im- 



portance of lattice distortions ior the stability of anti- 
ferromagnetism, Feiner and OlesE-3 favoured a purely elec- 
tronic mechanism. 

Our calculation points out that both parameters, U /Jh 
and g, can drive a transition from ferromagnetic (FM) to 
anti-ferromagnetic (AFM) order. Figure 0(c) shows the 
phase diagram of the electronic model without electron- 
phonon interaction, i.e., g — 0. We assume t = 0.4 eV 
and t/tn = 3 for the hopping integralstil and characterise 
the magnetic phases according to the total spin S tot of 
the ground-state of the four site cluster. Parameters in 
the range U < 5Jh implicate that for two neighbouring 
sites a d 5 d 3 configuration becomes the ground-state in 
favour of the d 4 ( 5 E)d 4: ( 5 E) configuration. This is incom- 
patible with the situation in the manganites and conse- 
quently the electronic Hamiltonian H e \ is not applicable 
for these values of U and J/,. Starting from the FM 
phase with 5J/j < U < 9.2Jh both, increasing U or g, 
change the magnetic order of the ground-state to AFM 
[Figure ||(a) and (b)]. The magnetic transition is accom- 
panied by a change in the corresponding orbital order, 
which we identify by the local expectation value (ng — n £ ) 
and the orbital exchange correlation (<Ti<Ti + $) between 
neighbouring sites. Here the Pauli matrices (af)^ oper- 
ate on the orbital degree of freedom, \i, v € {0, e}, at the 
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FIG. 6: Orbital-exchange [panels (a) and (b)] and spin-orbital 
[panels (c) and (d)] correlations between nearest and next- 
nearest neighbouring sites at variable electron-phonon cou- 
pling g and Coulomb interaction U. 



site i. The panels (a) and (b) of the Figures || and || il- 
lustrate the transition from staggered to uniform orbital 
order. 

In view of the distinct driving interactions (U or 
g) both transitions appear to be very similar. How- 
ever, we observe a significant difference, if we study the 
(de)coupling of spin and orbital degrees of freedom. The 
latter has been a rather controversial issue in the case of 
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the Kugel-Khomskii modelt 



U = 6.0 eV, J„ = 0.7 eV, t = 0.4 eV, t/t =3, o> = 70 meV 



KSiSiXTT + hi-rf + h 



+ tt*-o)(tf-o)-l 



(13) 



which contains the same kind of spin-orbital interac- 
tions SiSjT^Tj as our Hamiltonian. Using equation 
of motion approaches for the approximatp-,solution of 
the model (F|) Khaliullin and OudovenkcEZI decoupled 
spin and or jital degrees of freedom, while Feiner, Oles 
and ZaanenEa emphasised the role of mixed spin-and- 
orbital excitations. A correlation function, which makes 
it possible to distinguish between these two decoupling 
schemes, is the spin-orbital fluctuation (SiSi + ^ rfrf.g) — 
(SiSj+a) (rfrfig). The data we calculated for the com- 
plete manganite model ( fl^ ) is given in Figure |^(c) 
and (d). It indicates that this correlation is smaller by 
a factor of 3 to 5, if phonons are responsible for the 
FM to AFM transition (open symbols). On the other 
hand, starting within the AFM phase increasing electron- 
phonon interaction clearly suppresses the spin-orbital- 
fluctuation [Figure ||(c), bold symbols]. This behaviour 
is, of course, crucial for effective theories that are based 
on such decoupling schemes. 

To characterise the behaviour of the lattice we calcu- 
late the correlations (qi^ a Qj,a) (ot G {9,e}) between the 
elongations of the Jahn- Teller modes qg and q e at neigh- 
bouring sites i, j. The bond lengths 



(14) 



define additional significant quantities. In particular we 
considered the expectation values (q^s) and (qf s ) — (qi t s) 2 
with 5 G {x,y, z}. Figure |?] shows that within the FM 
domain the lattice is undistorted but the q e vibrations are 
correlated along the (tv, tt) direction in reciprocal space. 
At intermediate values of g/ui f=a 1 a finite x-y-distortion 
develops, which is also reflected in spin and orbital corre- 
lations [e.g., the spin-orbital fluctuation in Figure |](c)]. 
Such type of asymmetry seems surprising in view of the 
invariance of the cluster with respect to diagonal reflec- 
tions (dot-dot-dashed in Figure ||). However, since this 
symmetry is not taken into account in the calculation the 
system is trapped in one particular linear combination of 
two degenerate ground-states and different expectation 
values in x and y direction can evolve. At larger electron- 
phonon coupling g only a finite qg distortion remains and 
the q £ modes are uncorrelated. 

In Figure || the change of orbital, spin and phonon cor- 
relations with increasing g is shown schematically. This, 
however, is a rather suggestive picture, which can not 
cover all the details of the various correlations. To select 
the appropriate graphical representation for the orbital 
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FIG. 7: Behaviour of the lattice with increasing electron- 
phonon coupling g: (a) expectation value, and (b) fluctua- 
tion of the bond lengths qs; (c) and (d) correlations of the 
Jahn- Teller modes qe and q E between different sites. 




FIG. 8: Evolution of lattice, spin and orbital correlations with 
increasing electron-phonon coupling g at doping x = 0. 



arrangement we study the E ® S-eigenstates of reduced 
orbital density matrices. On a bond (ij) these states are 
either anti-symmetric, 

\a) fj « flfy ® \e)j - \e)i ® \e) 3 ) , (15) 

or symmetric, 

\s(ip, cx (\<p)i ® \^) 3 + <g> \<p)j) , (16) 

with respect to the permutation of two sites. The an- 
gles <p, ip G C parameterise two generalised orbital states 
of the form \ip) = cos(ip)\9) + sin(ip)\e). At small g the 
staggered orbital order implies that the anti-symmetric 
orbital state is the most probable for each bond, i.e., it 
belongs to the largest eigenvalue of the density matrix. 
Since an anti-symmetric combination of two arbitrary 
generalised states \ip) and is always proportional to 
\a)ij the representation is undefined. Therefore we took 
into account also the second eigenstate of the density ma- 
trix, which is slightly less probable but necessarily sym- 
metric. The associated angles tp and ip define the orbital 
pattern shown in the left hand panel of Figure [§. For 
larger g the uniform orbital order allows for a calculation 
of the depicted orbital pattern directly from the most 
probable eigenstate of the density matrix. The spin ar- 
rangements are chosen such as to reflect the expectation 
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FIG. 9: Doping x — 0.25. (a) Total spin and kinetic energy vs. 
g; (b) Generalised orbital states \ip) surrounding empty sites; 
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increasing g 




FIG. 10: Evolution of lattice, spin and orbital correla- 
tions with increasing electron-phonon coupling g at doping 
x = 0.25. 



values of the Heisenberg interaction, (S;Sj), for nearest 
neighbour bonds (ij). Of course, the calculated ordering 
patterns can not be directly compared to experimental 
observations for the three dimensional compounds. Nev- 
ertheless, the orbital arrangement within the FM spin 
background closely resembles the data measured within 
the FM planes of the A-type AFM structured. 



B. Doping x = 0.25 

In view of the CMR effect the domain of low to inter- 
mediate doping (0.15 < x < 0.5) is certainly the most 
interesting one. Here ferromagnetism is stabilised by 
the double exchange interaction, which depends on the 
mobility of the charge carriers. However, if too strong 
electron-phonon coupling causes localisation of the holes, 
the spin order breaks down. For x — 0.25 our calculation 
of the ground-state properties clearly illustrates this co- 
incidence. Figure ^ (a) shows the total spin S tot of the 
cluster together with the expectation value of the kinetic 



kin 



Rs \ a i,° a . 



3,o- U i,6 



H.C. 



2 shows a kink near the 
This behaviour closely 



(17) 

Obviously, the change from FM to AFM spin correla- 
tions with increasing electron-phonon interaction g is di- 
rectly related to the charge mobility. In addition the 
transition is accompanied by the appearance of a lattice 
distortion in the x-y-plane [see Figure |(c)], which is in- 
teresting also in view of the lattice dynamics. Whereas 
one of the elongations, {q v ), grows linearly in g, the as- 
sociated fluctuation (q 2 ) — (q y ) 
FM to AFM transition point. 

reminds the data collected by Booth et al.c2l using x 
ray-absorption fine-structure measurements (XAFS). For 
Lai^zCa^MnOa these authors observed a similar rise in 
the Mn-0 bond length variance a 2 near the critical tem- 
perature of the transition from the ferromagnetic metallic 
to the paramagnetic insulating phase. 

Orbital correlations are again extracted from the eigen- 
states of the reduced density matrices for a single bond. 
We concentrate on the environment of a hole, which is 
characterised by the E&A2 eigenstates, i.e., by states de- 
scribing a Mn 3+ site neighbouring a Mn 4+ site. Each of 
these states can be understood as the product \(p)i <g> \a,2) j 
of a generalised orbital state \ip) and the basis state \a2) 
[see Eq. @], In Figure ^|(b) the corresponding angle 
ip of the most probable eigenstate of the density ma- 
trix is given for bonds in x and y direction. The data 
for weak and strong electron-phonon coupling g is trans- 
lated into the orbital arrangement presented in Figure [lC|. 
Increasing g reduces nearest-neighbour correlations and 
forces the orbital order to depend locally on the dorak 
nant electron-phonon interaction. The orbital polaroncll 
pattern disappears, if the charge carrier is trapped com- 
pletely. 



C. Doping x = 0.5 

At doping x = 0.5 the picture is more involved. Again 
we consider model parameters, which yield ferromagnetic 
spin correlations for vanishing electron-phonon coupling, 
g = 0. The strong Coulomb interaction U then leads to 
the formation of a charge density wave with ordering vec- 
tor (7r, 7r) in reciprocal space. Increasing electron-phonon 
coupling g amplifies this trend, and at large g the trap- 
ping of the carriers gives rise to charge order. In Fig- 
ure |ll| (c) this is illustrated with the density-density cor- 
relation between different lattice sites i and j. From the 
model Hamiltonian (||) it is obvious that if charges tend 
to maximise their mutual distance the particular anti- 
ferromagnetic component of the Heisenberg spin interac- 
tion, which is proportional to t 2 /Jh, gains importance. 
Consequently FM order is unstable at much lower g, if 
compared to the case x = 0.25. The FM to AFM tran- 
sition is not connected to charge localisation and causes 
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FIG. 11: Doping x = 0.5. (a) Total spin and kinetic energy vs. 
g; (b) Complex orbital states on diagonal bonds; (c) Density- 
density correlations, (d) Lattice distortion (qg) in x and y 
direction. 




FIG. 12: Evolution of lattice, spin and orbital correlations 
with increasing electron-phonon coupling g at doping x = 0.5. 



only a tiny jump of the kinetic energy [see Fig. [Tl] (a)]. 

The calculated orbital correlations are interesting as 
well. On a diagonal bond the most probable E ® E- 
eigenstate of the reduced orbital density matrix is sym- 
metric for all considered values of g, \s(ip,ip)}. Note 
however, that the associated angles are complex conju- 
gated ip = ip* S C within both the FM and parts of the 
AFM phase [see Figure |Tl|(b)]. Here exa&Ldlagpnalisa- 
tion supports recent mean-field studiescjaOEJ, which 
obtain different patterns of complex orbital ordering, in 
particular within the ferromagnetic conducting phase. At 
large g charge localisation promotes real orbital states. 

In Figure |l2| we summarise the behaviour of the lattice, 
the spins and the orbitals. With growing g the cluster- 
expands isotropically until a finite x-y-distortion devel- 
ops at the FM to AFM transition [see also Figure pT|(d)]. 



Further increase of g causes the carriers to localise and 
the distortion to disappear. The spin background un- 
dergoes a change from FM order to different types of 
AFM order. Of course, within the four site system we 
are not able to detect more complicated spin arrange- 
ments like the CE-type order observed experimentally^. 
However, the sensitivity of the system to small changes in 
the model parameters is clearly visible. Note that the de- 
picted orbitals represent the amplitude of the underlying 
complex states \<p). 



IV. CONCLUSION 

In the present work we have derived a microscopic 
model for doped CMR manganites which includes the dy- 
namics of charge, spin, orbital, and lattice degrees of free- 
dom on a quantum mechanical level. Using exact diago- 
nalisation techniques we have studied how the electron- 
lattice interaction affects short range spin and orbital cor- 
relations. An observation, which is important for the un- 
derstanding of the undoped compounds, is the suppres- 
sion of the spin-orbital coupling with increasing electron- 
phonon interaction. For the weakly doped compounds we 
have demonstrated the direct relationship between the 
trapping of charge carriers and the breakdown of ferro- 
magnetism. In addition, we have shown that changes 
in the spin correlations are reflected in dynamic lattice 
correlations. At intermediate doping we find that the 
system depends on a subtle balance of double-exchange, 
superexchange and electron-lattice interaction. The lat- 
ter enhances the tendency for charge ordering which in 
turn affects spin and orbital order. Besides the calcula- 
tion proves that complex orbital states can be a suitable 
approximation for the description of orbital correlations. 
Although the calculated data is not quantitatively com- 
parable to real, three dimensional materials, the exact 
diagonalisation of even a small system provides detailed 
insight into correlations and driving interactions behind 
the rich phase diagram of the manganites. In addition, 
the exact results may support the development of ap- 
proximate theories. 

We thank L.F. Feiner, F. Gohmann, D. Ihle, G. Khal- 
iullin, and J. Loos for valuable discussions. The 
grant of computational resources by HLR Stuttgart and 
NIC Jiilich, as well as financial support by the Deutsche 
Forschungsgemeinschaft and the Czech Academy of Sci- 
ences is acknowledged. 



1 R. M. Kusters, J. Singleton, D. A. Keen, R. McGreevy, 
and W. Hayes, Physica B 155, 362 (1989). 

2 R. von Helmolt, J. Wecker, B. Holzapfel, L. Schultz, and 
K. Samwer, Phys. Rev. Lett. 71, 2331 (1993). 

3 S. Jin, T. H. Tiefel, M. McCormack, R. A. Fastnach, 
R. Ramesh, and L. H. Chen, Science 264, 413 (1994). 



4 J. M. D. Coey, M. Viret, and S. von Molnar, Adv. Phys. 
48, 167 (1999). 

5 Y. Tokura and Y. Tomioka, J. Magn. Magn. Mater. 200, 
1 (1999). 

6 E. Dagotto, T. Hotta, and A. Moreo, Physics Reports 344, 
1 (2001). 



8 



7 C. Zener, Phys. Rev. 82, 403 (1951). 

8 P. W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 
(1955). 

9 K. Kubo and N. Ohata, J. Phys. Soc. Jpn. 33, 21 (1972). 

10 A. J. Millis, Nature 392, 147 (1998). 

11 S. Ishihara, J. Inoue, and S. Maekawa, Physica C 263, 130 
(1996). 

12 R. Shiina, T. Nishitani, and H. Shiba, J. Phys. Soc. Jpn. 
66, 3159 (1997). 

13 R. Maezono, S. Ishihara, and N. Nagaosa, Phys. Rev. B 
57, R13993 (1998). 

14 L. F. Feiner and A. M. Oles, Phys. Rev. B 59, 3295 (1999). 

15 J. S. Griffith, The Theory of Transition-Metal Ions (Cam- 
bridge University Press, Cambridge, 1971). 

16 P. W. Anderson, Phys. Rev. 115, 2 (1959). 

17 K. I. Kugel and D. I. Khomskii, JETP Lett. 15, 446 (1972). 

18 A. Weifie, J. Loos, and H. Fehske, Phys. Rev. B 64, 054406 
(2001). 

19 Y. Tanabe and S. Sugano, J. Phys. Soc. Jpn. 9, 766 (1954). 

20 J. Zaanen and G. A. Sawatzky, J. Solid State Chem. 88, 8 
(1990). 

21 A. E. Bocquet, T. Mizokawa, T. Saitoh, H. Namatame, 
and A. Fujimori, Phys. Rev. B 46, 3771 (1992). 

22 T. Mizokawa and A. Fujimori, Phys. Rev. B 51, 12880 
(1995). 

23 M. Quijada, J. Cerne, J. R. Simpson, H. D. Drew, K. H. 
Ahn, A. J. Millis, R. Shreekala, R. Ramesh, M. Rajeswari, 
and T. Venkatesan, Phys. Rev. B 58, 16093 (1998). 

24 Y. E. Perlin and M. Wagner, eds., The Dynamical Jahn- 
Teller Effect in Localized Systems, no. 7 in Modern Prob- 
lems in Condensed Matter Sciences (North-Holland, Ams- 
t erdam, 1984). 



:-!() 



T |. Holstein, Ann. Phys. (NY.) 8, 325 (1959). 



29 J. K. Cullum and R. A. Willoughby, Lanczos Algorithms 
for Large Symmetric Eigenvalue Computations, vol. I & II 
(Birkhauser, Boston, 1985). 

E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 
(1955). 

31 F. Moussa, M. Hennion, J. Rodriguez-Carvajal, H. Moud- 
den, L. Pinsard, and A. Revcolevschi, Phys. Rev. B 54, 
15149 (1996). 

32 D. D. Sarma, N. Shanthi, S. R. Barman, N. Hamada, 
H. Sawada, and K. Terakura, Phys. Rev. Lett. 75, 1126 

(1995) . 

33 W. E. Pickett and D. J. Singh, Phys. Rev. B 53, 1146 

(1996) . 

34 S. Satpathy, Z. S. Popovic, and F. R. Vukajlovic, Phys. 
Rev. Lett. 76, 960 (1996). 

35 I. Solovyev, N. Hamada, and K. Terakura, Phys. Rev. Lett. 
76, 4825 (1996). 

36 K. I. Kugel and D. I. Khomskii, Sov. Phys. JETP 37, 725 
(1973). 

37 G. Khaliullin and V. Oudovenko, Phys. Rev. B 56, R14243 

(1997) . 

38 L. F. Feiner, A. M. Oles, and J. Zaanen, J. Phys. Condens. 
Matter 10, L555 (1998). 

39 Y. Murakami, J. P. Hill, D. Gibbs, M. Blume, I. Koyama, 
T. Tanaka, H. Kawata, T. Arima, Y. Tokura, K. Hirota, 
et al, Phys. Rev. Lett. 81, 582 (1998). 

40 C. H. Booth, F. Bridges, G. H. Kwei, J. M. Lawrence, 
A. L. Cornelius, and J. J. Neumeier, Phys. Rev. Lett. 80, 
853 (1998). 

41 R. Kilian and G. Khaliullin, Phys. Rev. B 60, 13458 
(1999). 

D. Khomskii (20 00), URL |http : //arXiv . org/ab"s/| 



M. N. Iliev, M. V. Abrashev, H.-G. Lee, V. N. Popov, Y. Y. 
Sun, C. Thomsen, R. L. Meng, and C. W. Chu, J. Phys. 
Chem. Solids 59, 1982 (1998). 

27 W. Reichardt and M. Braden, Physica B 263-264, 416 
(1999). 

28 A. Weifie, H. Fehske, G. Wellein, and A. R. Bishop, Phys. 
Rev. B 62, R747 (2000). 



cond-mat/0004034 . 

4o A. Takahashi and H. Shiba, J. Phys. Soc. Jpn. 69, 3328 
(2000). 

44 R. Maezono and N. Nagaosa, Phys. Rev. B 62, 11576 
(2000). 

45 J. van den Brink and D. Khomskii, Phys. Rev. B 63, 



140416 (2001). 



